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Abstract 

Recent gyrokinetic stability calculations have revealed that the spherical tokamak is susceptible to tearing 
parity instabilities with length scales of a few ion Larmor radii perpendicular to the magnetic field lines. Here 
we investigate this 'micro-tearing' mode in greater detail to uncover its key characteristics, and compare it with 
existing theoretical models of the phenomenon. This has been accomplished using a full numerical solution of the 
linear gyrokinetic-Maxwell equations. Importantly, the instability is found to be driven by the free energy in the 
electron temperature gradient as described in the literature. However, our calculations suggest it is not substantially 
affected by either of the destabilising mechanisms proposed in previous theoretical models. Instead the instability 
is destabilised by interactions with magnetic drifts, and the electrostatic potential. Further calculations reveal that 
the mode is not significantly destabilised by the flux surface shaping or the large trapped particle fraction present in 
the spherical tokamak. Its prevalence in spherical tokamak plasmas is primarily due to the higher value of plasma 
/3, and the enhanced magnetic drifts due to the smaller radius of curvature. 
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O ■ 1 Introduction 

■ 

£>-y The linear tearing instability, described by Furth, Killccn, and Rosenbluth [T|, progresses by relaxing a sheared 
' magnetic field into a lower energy magnetic island structure. Importantly, the energy released from this process 
l Ph. decreases as a function of mode number, which prevents the formation of shorter wavelength tearing instabilities. 

\ However, later calculations by Hazeltine et al [2], demonstrated that the linear tearing instability could also be driven 
t-H ■ by an electron temperature gradient in slab geometry. This new mechanism may enable the tearing instability to exist 
^ \ at higher mode numbers, giving rise to the term 'micro-tearing mode.' 

For a given toroidal mode number n, the distance between adjacent rational surfaces is Sr = where q' is the 
radial derivative of the safety factor q. Consequently, the chains of islands formed by micro-tearing instabilities are in 
■ close proximity, and may overlap producing a stochastic field line structure. When this occurs, the electrons are no 
( ^-^ longer confined to equilibrium flux surfaces, and are able to move radially by moving along the perturbed field lines. 

I ■ This process can produce significant electron heat transport, which could explain the high ratios of Ti/T e measured 
y—i \ in a number of spherical tokamak (ST) plasmas (including those from START and NSTX [3111]). 

The drive mechanism for this micro-tearing instability can be explained in terms of the parallel thermal force [5]- 
This force is due to the different frictional forces experienced by electrons travelling in opposite directions along a 
V~j ■ parallel temperature gradient, and has a magnitude given by, 
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where n e and T e are the electron density and temperature respectively. The parallel current resulting from this force 
will produce a magnetic field perturbation B. Instability then arises if this field perturbation aligns with an equilibrium 
electron temperature gradient VT e . The later paper by Hassam [5], uses a second order Chapman-Enskog expansion 
of the fluid equations to recover a similar solution. This work notes that the time dependent thermal force is also 
important to the drive mechanism of the micro-tearing instability. (This time-dependent force is not present in first 
order Chapman-Enskog fluid theories, so the Braginskii equations [B] are unable to recover the instability). 

Drake and Lee [7] expand on the work of Hazeltine et al [5] , by carefully defining three collisional limits in which 
the slab/cylindrical micro-tearing instability can be analytically treated: the collisional regime, where v e » w; the 
semi-collisional regime, defined by D e > ui and v e > k'^v't h /u) (nb. perturbations take the form £ = £ e n me -»<p-< t, < 
such that fc|| = (m — nq)/Rq)\ and finally the collisionlcss regime where v e « lu. In these expressions, Vth c is the 
thermal velocity of electrons, uj is the mode frequency, k\\ is the wavenumber of the mode parallel to the equilibrium 
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field lines, and v e = ^^""3/2 '"1/2 determines the collision rate experienced by electrons. Hot tokamak plasmas typically 
have \> e < lo* (where w* = —k y T e /(eBn)dn/dr is the electron diamagnetic frequency), which best agrees with the 
collisionless limit of Drake and Lee. Importantly the thermal force vanishes in this limit, meaning the collisionless 
micro-tearing mode is stable in slab geometry [7J|S] suggesting the micro-tearing mode may not be important in fusion 
plasmas such as ITER. 

However, the paper by Catto and Rosenbluth [9] demonstrates another drive mechanism for the tearing instability 
in the ITER relevant banana regime {D e < io\, e where ujbe is the electron bounce frequency). In this limit electrons 
close to the trapped-passing boundary can easily scatter between the trapped and passing populations. This increases 
the effective collision rate in a layer around the trapped-passing boundary and allows a destabilising current to flow. 
Catto and Rosenbluth [5] have shown this can destabilise the tearing mode in the regime v e jt < lu. In similar spirit, 
Connor et al [10] show the trapped particle driven mode is stable when D e /e > 10 > D e . 

Importantly, gyrokinetic simulations have identified the micro-tearing mode in a number of ST H-mode plasmas 
[HJ EH 03] ■ Here we build on these previous calculations with a closer investigation of the instability. In particular, the 
key characteristics of the instability are uncovered, and compared with existing theoretical models of the phenomenon. 
Section [2] of this paper describes the equilibrium model we have considered, and also the system of equations employed 
in the study of the micro-tearing mode. In Section [3] the key features of the instability are discussed, including the 
sensitivity to temperature/density gradients, beta, collisionality, geometry, trapped particles, and also the importance 
of the electrostatic potential and the ion response. These results confirm the micro-tearing modes found in MAST are 
not well described by the existing literature, and reveal the processes that are crucial for instability. Our conclusions 
are given in Section 01 



2 Plasma Model 

The plasma considered is MAST discharge #6252, which is a well diagnosed Elmy H-modc plasma with neutral beam 
heating. An equilibrium reconstruction for this plasma has been calculated during a steady level of the plasma current, 
and just before the neutral beams were switched off. (The key plasma parameters are listed in Table [TJ and a detailed 
discussion of the equilibrium calculation has been given in a previous publication In particular we study the 

ip n = 0.4 poloidal flux surface of this equilibrium (ip n = ^ where ip is the poloidal flux), since this was found 

to be highly susceptible to the micro-tearing mode. (The important plasma parameters for this surface are given in 
Table [5]) ■ A significant difference between the plasma model employed here and that used in previous work , is 
the reduction in the number of kinetically modelled species from five to just two (main ions and electrons). This has 
significantly improved the tractability of our calculations without significantly affecting the results (see Figure [TJ). 
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Table 1: Global parameters for the MAST equilibrium at t = 0.265s in MAST discharge #6252. Here B is the 
toroidal magnetic field on the magnetic axis, I p is the plasma current, R and a are respectively the major and minor 
radius of the last closed flux surface, R m ag is the major radius of the magnetic axis. Pbeam is the applied neutral 
beam heating, and finally Z e ff = Q 2 s n s /n e (where s labels the ion species and Q s is the electric charge of the ion 
species.). 
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Table 2: Important parameters for the poloidal flux surface ip n = 0.4. T e and Tj are the electron and ion temperatures, 

P is the plasma pressure, while n gives the number density for both ions and electrons. The normalised gradient of 

, \ 1/2 , \ 3/2 

a quantity x is defined as a/L x = -j^, while s = The parameters i>i — D e Q 2 Z e fj f^r-J (tM > an< ^ 

D e = (^"3/2 '"i/a ' determine the collision rates of ions and electrons (respectively) and arc both normalised to v t \ li / a. 
The quantity r is the half diameter of the flux surface across its midplane, while pi is the Larmor radius of thermal 
ions. 
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Figure 1: Growth rate (7) and real frequency (w) of micro-tearing modes on the ip n = 0.4 surface (both are normalised 
to Vthi/a). These arc plotted as functions of the normalised wave number k y pt (see Appendix A for the definition of 
the x and y coordinates used in GS2). The solid line represents the two species model, while the dashed line represents 
the five species model. 



The linear gyrokinetic-Maxwell equations [151 116] have been employed to study the micro-tearing instability present 
in this equilibrium. These equations provide a kinetic description for the plasma in the limit, 
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where fl is the gyro-frequency, w is the oscillation frequency of the field perturbations, k± and fen are the perpendicular 
and parallel wavenumbers, and L is the scale length of variations in the equilibrium distribution function. The linear 
gyrokinetic equation for a Fourier mode kj_ is given by the expression below. 
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where R = r— is the guiding centre coordinate; the nonadiabatic distribution function is given by g s = f a ~, j, — 

(where Q s is the particle charge, <j> is the electrostatic potential, f s is the perturbed distribution function, and Fq 



is the Maxwell-Boltzmann distribution function); finally £kj_ = <Pk±Jo{ z ) — v \\A\\ t k ± Jo(z)+ B^^ ± Ji{z) 
z = kj_v±/Q. Particle collisions are modelled using a Lorentz collision operator, 
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where v is the particle collision rate, £ = v\\jv is the pitch angle variable, and |r denotes a partial derivative at constant 
R. This operator calculates the pitch angle scattering of particles due to collisions, which enables an accurate model 
of trapped particle interactions. It also conserves energy and particle numbers. However it does not model the energy 
scattering involved in collisions, and does not include momentum conservation. The diffusive term proportional to 
k\ is found to have little influence on the linear micro-tearing mode [12] . However, preliminary work suggests it 
may have a more important effect in nonlinear simulations of the instability [17j . Both ion-electron collisions and 
electron-electron collisions are included in the electron collision rate Veiy), while only ion- ion collisions are included 
in the ion collision rate Vi{v). However, it is found that electron- ion collisions are by far the most important of these 
collisional processes. 

The gyrokinetic-Maxwell equations have been solved numerically using the initial value code GS2 [T5]. This code 
assumes the equilibrium temperature and density scale lengths are constant across the radial simulation domain, and 
since the radial simulation domain is generally quite small for the micro-tearing instability in MAST, this approx- 
imation is reasonable (see reference Two further assumptions have also been made to improve the numerical 
tractability of these calculations . First the parallel magnetic perturbations (B|| = V x Aj_) have been neglected, 
which, despite the high plasma /3, have little effect on the micro-tearing mode (see Section [575]) . Second, the plasma 
has been approximated as quasi-neutral. (Results obtained using the quasi-neutrality condition were found to be in 
good agreement with those obtained from Poisson's equation.) 



3 Results 

Figure [2] is a plot of the parallel electric field E\\ , the electrostatic potential <fr, and the parallel vector potential A\\ 
resulting from the micro-tearing instability at k v pi = 0.5. Each quantity is plotted as a function of the ballooning 
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angle 9 (which labels the distance along a magnetic field line using the poloidal angle). Importantly, the electrostatic 
potential and parallel electric field are highly extended along the field lines due to the rapid parallel motion of electrons. 
The parallel vector potential is confined to a much smaller range in 9, which, in real space, produces a radially extended 
function. This corresponds to the "constant -0" approximation of linear tearing theory pQ . 
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Figure 2: Eigenfunctions for the electrostatic potential (f>, the parallel vector potential A\\, and the parallel electric 
field £j| plotted as functions of the ballooning angle 9. The real part of each eigenfunction is represented by a solid 
line, whilst the imaginary parts are represented by dotted lines. 

The tearing nature of this instability becomes apparent when Poincare sections are plotted of the perturbed 
magnetic field. In GS2 this is achieved using a flux tube simulation domain (see reference [19] for details), which is a 
long thin tube aligned parallel to the magnetic field lines. The parallel length extends from 9 = [— tt, it] (as opposed to 
the arbitrarily large domain of the ballooning representation), while the perpendicular dimensions are a few Larmor 
radii. Parallel vector potential data from GS2 is then used to determine the path of selected magnetic field lines 
through the flux tube. Once these paths have been established, a perpendicular slice is taken through the tube at 
constant poloidal angle 9, and the points where the field lines intersect with this slice are plotted. 
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Figure 3: Poincare plots for the micro-tearing instability with k y pi=0.66 at 9 = 0, where the x and y axis are given in 
units of pi. The amplitude of the perturbation increases from left to right, so that the middle plot has a magnitude 
five times greater than the left hand plot, and four times smaller than the right hand plot. The right hand plot also 
includes a contour diagram of the parallel current J\\ . 

The Poincare plots in Figure |3] were produced using A\\ data for a single micro-tearing mode with wavenumbcr 
kyPi — 0.66. Each plot is given on an (x,y) grid, where x and y are the perpendicular coordinates in the flux tube, 
and are both normalised to the ion Larmor radius pi. (The definitions of x and y are given in Appendix A). A second 
horizontal axis labels the equivalent q values for the x axis. Importantly, the field lines are tracked from the same 
starting points in each plot, while the amplitude of the magnetic field perturbation (due to the instability) increases 
from left to right. A magnetic island structure appears at the rational surface q = || of each plot which increases 
in size and complexity as the amplitude of the perturbation increases. The current plot in Figure [3] reveals the 'O' 
point of the island is at the current minimum, whilst the 'X' point of the separatrix is at the current maximum. This 
current perturbation occurs in a narrow region around the resonant surface, and corresponds to the current layer of 
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linear tearing mode theory [TJ. 



3.1 Drift-wave Characteristics of the Micro-tearing Mode 

The thermal force and trapped particle drive mechanisms for the micro-tearing instability both require a finite electron 
temperature gradient from which to draw free energy. The left hand side of Figure 0] is a plot of growth rate as a 
function of a/LTe for the MAST micro-tearing instability with k y pi = 0.5. (Note, the ion temperature gradient 
was simultaneously altered to maintain a constant pressure gradient.) The instability is clearly dependent on the 
presence of a finite electron temperature gradient, as predicted by Hazcltinc et al [2]. Further calculations have shown 
the growth rate also varies with the electron density gradient, although a finite density gradient is not vital to the 
instability. 
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Figure 4: Growth rates and real frequency of micro-tearing modes as a function of the inverse electron temperature 
scale length a/Lx e - For comparison, the dashed line in the real frequency plot represents the analytic prediction 
w*[l+T) e /2]. Note that the inverse electron density scale length was kept constant, so r\ e — varies proportionally 



with a/L 



Previous theories of the micro-tearing mode often predict a linear relation between the equilibrium electron tem- 
perature and density gradients and the real frequency of the mode. For example the paper by Catto and Rosenbluth 
[9] predicts ui ~ + r] e /2] (where r\ e = jf^), while numerical calculations by Gladd et al [20j produce real frequen- 



ce] . Both of these expressions have a qualitative 



cies which roughly fit the expression, uj ~ w*[l + 1 + 0.6 In {j^j 

agreement with calculations for the MAST instability (see Figure H] for the dependence of the real frequency on the 
electron temperature gradient). However, the frequency of the MAST instability also depends on the collision rate, 
which is closer to the predictions of Gladd et al. 

Curiously there is a minimum in the growth rate of Figure 0] at a/LTe ~ 4.5, which coincides with a discontinuity 
in the frequency plot. This signifies a mode switching between two different micro-tearing instabilities, one which is 
unstable for < a/LTe < 4.5 (which we will call mode A), and another which is unstable for a/LTe > 4.5 (which we 
will call mode B). It is not obvious why mode A should be stabilised above a/LTe ~ 4.5 since the temperature gradient 
provides the free energy for the instability. A plausible explanation could be the dependence of the real frequency on 
the electron temperature gradient. If a resonance is important (for example with the magnetic drift frequency), then 
increasing \uj\ may disturb this resonance, and cause mode A to be damped. Interestingly, \uj\ sharply decreases across 
a/LTe ~ 4.5. This means mode B has a real frequency similar to the real frequency of mode A at its most unstable, 
and supports the idea that a lower \uj\ produces a destabilising resonance in the plasma. 

3.2 The Role of Collisions 

It was discussed earlier that Drake and Lee [7] have defined three specific categories of tearing mode depending on 
the level of collisionality. However, it is not obvious which of these three categories the MAST instability belongs to. 
Consider Figure [SJ which shows the growth rate of the MAST micro-tearing instability as a function of collisionality. 
The instability is not collisionless, since the growth rate is clearly dependent on v e , and instability requires finite v e 
(Figure [5]). It is also not collisional or semi-collisional, since the collision rate is too low: 0.05 < v e /u) < 2. In fact 
the instability seems to lie somewhere between the collisionless and semi-collisional regimes. Importantly, the MAST 
instability is unstable for v e ~ w. This regime has been addressed by Gladd et al [2D] and D'Ippolitto et al [3T] using 
a modified Krook collision operator [33J[53]. Conversely, at lower collisionalities, the MAST instability is closer to the 
limit treated in the trapped particle calculations of Catto and Rosenbluth [9], where u> > u e /e. 

3.2.1 The Energy Dependence of the Collision Operator 

Importantly, the thermal force drive mechanism of Hazcltinc et al [2], and the trapped particle drive mechanism 
of Catto and Rosenbluth [3], both require an energy dependent collision operator. Consequently we have tested 
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whether the MAST instability is sensitive to the energy dependence of collisions. This test involved replacing the 
velocity /energy dependent collision rate of the electrons with the collision rate of electrons moving at their thermal 
velocity, 

u e (v) ->• v e (v t h e ) (5) 

Figure [5] is a plot of growth rate and real frequency versus the electron collisionality P e , using both the energy 
dependent collision operator and the energy independent collision operator. It is clear the instability is not strongly 
affected by the change in collision operator, and thus cannot be strongly dependent on either of the drive mechanisms 
mentioned in the literature. 

However, in the absence of magnetic drifts (see Section 13. 5|) energy dependent collisions become more important 
to the MAST instability. This makes contact with the numerical calculations of Gladd et al [50] and D'Ippolitto et al 
j21j . who find a micro-tearing instability in slab geometry, which is reliant on energy dependent collisions, and occurs 
in the regime v e ~ uj. 
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Figure 5: Plots of growth rate and real frequency for micro-tearing instabilities as a function of electron collisionality 
D e , for two different collision operators. The solid lines represent the energy dependent collision operator, while the 
dashed lines represent the energy independent collision operator from equation (J5J). 
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3.3 The Role of the Ion Response 

The appropriate form for the ion response in tearing mode calculations has been widely discussed in the literature. 
Early calculations by authors such as Drake and Lee [7J, utilise magnetised or unmagnetised approximations for the 
ion response, depending on the predicted width of the current layer d (a magnetised response is used if d ^> p, while 
an unmagnetised response is used for d <C Pi). Later calculations by Cowley et al [8], implimented the effect of finite 
ion Larmor radius (FLR) effects more completctly, in the limit d < pi. To discern the effect ion physics has on the 
MAST micro-tearing instability, a comparison has been made between calculations using the full ion response, and 
calculations using an unmagnetised ion response, 

ed> , . 

fH = ~7^F (6) 

where hi is the perturbed ion density. The results of these calculations are presented in Figure |5J where the growth 
rate and real frequency of the instability (at k y pi = 0.5) are plotted against collision frequency. It is clear the full 
ion response has a limited influence on the instability, contrasting with the calculations of Cowley et al [5], who find 
FLR effects have a strong stabilising effect on the collisionless and semi-collisional tearing instabilities. However, 
this apparent contradiction can be explained by considering the current layer produced by the MAST micro-tearing 
instability. Figure [7] is a contour plot of the current perturbation produced by the k y pi = 0.5 micro-tearing mode at 
6 = 0. The current perturbation is concentrated in a small layer around x = 0, and here we define the width of this 
layer, d, using the following expression, 

Loo \ J \\\ dxd v 

Figure [7] reveals d ~ pi throughout the unstable collisional range, suggesting the approximation d < pi, made by 
Cowley et al [8], is not applicable to the MAST instability. (In the limit d < pi the current layer itself is largely 
unmagnetised, but the full ion response is still important outside the current layer. Indeed it is finite Larmor radius 
(FLR) effects outside of the current layer which Cowley et al find has a stabilising effect on the instability. On the 
other hand the MAST instability has d/pi ~ 1 so ion FLR effects ought to be important both inside and outside of 
the current layer.) 
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Figure 6: Plots of growth rate and real frequency versus collisionality for calculations with a full ion response (solid 
lines), and a Boltzmann response for the ions (dashed lines). 




Figure 7: (Left) Cross section of the parallel current Jn in the flux tube at 8 = 0, and at k y pi = 0.5. The dashed 
lines show the width d of the current layer which is defined in equation ([7]). (Right) The current width d is plotted as 
a function of collisionality 9 e , using full ion physics (solid line), and also a Boltzmann response for the ions (dashed 
line) . 



3.4 Micro-tearing Instabilities in s — a Geometry 

While numerical equilibria provide a realistic reconstruction of the experimental situation, they are somewhat inflexible 
when one wishes to modify equilibrium parameters. Thus in order to study the effect that trapped particles, magnetic 
drifts and aspect ratio have on the micro-tearing instability, it is helpful to introduce an analytic equilibrium. Here 
we have constructed a shifted-circle, s — a equilibrium |24) where the flux surfaces are specified by, 



R = Rq (r) + rcos6 

Z = rsm6 (8) 

where r is the minor radius of the flux surface, R and Z are the horizontal and vertical coordinates that specify the 
flux surface, and Ro(r) is the major radius of the flux surface (the s — a model allows this to be a function of the flux 
surface label r). In this model the total magnetic drift velocity for electrons is given by, 

2«jj + "1 R x B 

v R + WB = - " 2Q R2B (9) 

Trapped particles can be included in the s — a model by using the following expression for the variation of the 
magnetic field strength, 



B = — ^— rrr (10) 

1 + ecos(6») v ' 

where e = r/ Rq is the inverse aspect ratio of the flux surface and controls the fraction of trapped particles. 

In Figure [8] the growth rates and real frequencies obtained using the s — a model are compared with those from 
the numerical MAST equilibrium. (This is only a rough comparison since the self consistency of the s — a model 
deteriorates at high (3 [25]). The same values of q, /?, s, q/Lt, a/L n , Rq and r/R are used in both models. It is 
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Figure 8: Growth rate and real frequency versus collisionality for four different plasma geometries. Calculations using 
the numerical reconstruction are given by solid lines. Results from the s — a model closest to this equilibrium (using 
r/Ro = 0.345 and Ro = 0.879) are given by dashed lines. The s — a calculations using r/Ro = 0.236 and Ro = 1.32 
are given by dotted lines. Finally s — a results using r/Ro = 0.118 and Rq = 2.64 are given by dash-dot lines. 

clear the micro-tearing mode is unstable in the s — a equilibrium, suggesting the flux surface shaping of the MAST 
equilibrium is not particularly important for instability. 

Using the s — a model, the major radius Ro and inverse aspect ratio e = r/Ro of a flux surface can be altered. In 
Figure[Hl the value of Ro has been changed keeping r fixed, and allowing r/Ro to vary. Importantly, all other quantities, 
such as the temperature and density scale lengths, are kept constant. This is roughly equivalent to changing the aspect 
ratio Ro/a of a tokamak, whilst staying on the same flux surface r. From Figure [5] it is clear that any increase in 
aspect ratio has a stabilising effect on the micro-tearing instability, although it is also apparent the mode remains 
significantly unstable if the major radius is modestly increased. This in turn suggests the instability might appear in 
moderate aspect ratio machines such as JET (although it will be shown later on that the instability also requires a 
large /?, which may prevent it from appearing at large aspect ratio). 

3.5 The Effect of Magnetic Drifts 

In the previous section an s — a model was utilised to vary Ro and e = r/Ro simultaneously. Consequently, both the 
magnetic drift velocity (which is inversely proportional to Ro) and the trapped particle fraction (which is given by 
e 1 / 2 ) were changed. In this section the s — a model is used to vary Ro while keeping e = r/Ro constant, allowing us 
to consider the effect of magnetic drifts in isolation. 

The role that magnetic drifts play in the micro-tearing instability has not been fully treated in the literature. The 
majority of papers use either slab or cylindrical geometry to help tractability [5J [7J [501 IB], while a smaller number 
attempt to understand the role that trapped particles have on the instability [HHUl], but ignore the effects that curvature 
and V-B drifts have on untrapped electrons. Figure [9] compares the growth rates and real frequencies for calculations 
with a number of different major radii, and with a range of collisionalities. These results show magnetic drifts have a 
significant destabilising effect on the instability, and it is found that the destabilisation comes predominantly via the 
untrapped electron species. 

It is tempting to suggest the instability is driven by bad curvature at the outboard midplane of the tokamak, which 
is the destabilisation mechanism for the toroidal ITG/ETG mode. However, in the case of the micro-tearing mode, the 
passing electrons must travel many poloidal orbits in an oscillation time of the mode (^), and would be very sensitive 
to the average magnetic drift rather than just the bad curvature at the outboard midplane. Also, calculations have 
been performed where a = — ^^f 3 has been increased in order to reduce the average curvature (this has been 
done inconsistently so that the driving gradients are kept constant). This was found to increase the growth rate of 
the micro-tearing instability, whereas bad curvature driven modes, such as the toroidal ITG mode, are expected to 
be stabilised. Thus, the destabilising influence of magnetic drifts appears to be due to another process (perhaps a 
resonance between the mode frequency and the magnetic drift frequency as suggested in Section |3~T|) . 

The results from Figure [0J show there is also an underlying slab drive mechanism for the micro-tearing instability 
in MAST, since the mode is unstable even when the magnetic drift velocity is zero. As we mentioned earlier, this 
additional mechanism appears to be closer to the slab calculations of Gladd et al [20], since the instability becomes 
much more sensitive to the energy dependence of the collision operator in the absence of magnetic drifts. 

3.6 Calculations without trapped particles 

It is also possible to investigate the role that trapped particles play in the micro-tearing instability using the s — a 
model. To do this we simply vary the ratio r/Ro whilst keeping Ro fixed. Figure [TU] compares the growth rates and real 
frequencies of calculations with a number of different values of r/Ro and a range of collisionalities. Importantly these 
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Figure 9: Plots of growth rate and real frequency versus collisionality using the s — a model with four different major 
radii at fixed inverse aspect ratio r/Ro- Solid lines represent Rq = 0.879, dashed lines represent Rq = 1.32, dotted 
lines represent Rq = 2.64, dash-dot lines represent Rq — oo (ie. cylindrical geometry). Note that k y pi = 0.5 for all of 
these calculations. 

collision frequencies are well below the electron bounce frequency which is w& = 12.6 (all frequencies are normalised 
to -^ J -) so the banana orbits of trapped electrons are important . 

Figure QJJ reveals that trapped particles are destabilising at lower collision frequencies, which is in general agreement 
with the predictions of Catto and Rosenbluth [3]. It is also apparent that trapped particles are most destabilising at 
smaller inverse aspect ratio, where the calculations of [9] are most accurate. 

Importantly the trapped particles have a stabilising effect at higher collision frequencies. This is probably because 
the trapped particles cannot respond to a parallel electric field, so in the absence of destabilising boundary layer 
effects, the trapped particles will significantly reduce the perturbed current in the plasma, and hence damp other drive 
mechanisms for the micro-tearing instability. Thus, we can see the stabilising effect becomes stronger as the trapped 
particle fraction increases, and is particularly strong for the aspect ratio of MAST. 




0.1 1 0.1 1 

Figure 10: Plots of growth rate and real frequency versus collisionality using the s — a model with four different inverse 
aspect ratios r/Ro at fixed major radius Rq. Solid lines represent r/Ro = 0.354, dashed lines represent r/Ro — 0.236, 
dotted lines represent r/R = 0.118, and dot-dash lines represent r/Ro = 0. 



3.7 The Electrostatic Potential 

Gladd et al [20| find the electrostatic potential has a strong destabilising influence on the micro-tearing instability when 
uj < v e . The importance of the electrostatic potential has also been tested for the MAST micro-tearing instability by 
completing calculations with and without the electrostatic potential for a wide range of collisionalities (see Figure ITT]) . 
It is immediately apparent that the electrostatic potential has a strong destabilising influence on the micro-tearing 
instability. Indeed if both the electrostatic potential and the magnetic drifts are turned off then the mode is stabilised. 

3.8 The Effect of (3 and B\\ 

Figure [T2l presents a scan of growth rate versus total (3. Two models of the plasma are compared, one including parallel 
magnetic perturbations (Bu) and one without. Importantly, these calculations keep (3' and VP constant while varying 
(3 (which is inconsistent), so the strength of magnetic perturbations in Ampere's Law is varied, while the equilibrium 
is kept the same. The micro-tearing instability clearly requires a significant j3 value for instability. Below /3 ~ 0.1 
(which is the (3 value of the MAST equilibrium) the growth rate decreases rapidly. This is understandable since the 
lower values of j3 make it more difficult for the magnetic field to be perturbed. However, above j3 <~ 0.15 the growth 
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Figure 11: Plots of growth rate and real frequency versus collisionality. Solid lines include <f>, whilst the dashed lines 
have = 0. 



rate decreases again, which is not understood at present. Finally, we note that B\\ is unimportant to the micro-tearing 
instability, even at /3 ~ 0.4. 



f- 0.05 




Figure 12: Plots of growth rate and real frequency versus plasma (3. Solid lines include parallel magnetic perturbations, 
dashed lines do not. 



4 Conclusion 

The tearing nature of the instability found in MAST has been confirmed by producing Poincare plots of the perturbed 
magnetic field. These reveal chains of magnetic islands on the mode rational surfaces, and it is around these rational 
surfaces that most of the current perturbation lies in a 'current layer.' The centre of each island coincides with current 
minima, while the 'X' points coincide with current maxima; and the width of each magnetic island grows as the 
amplitude of Au increases. These observations are in good agreement with linear tearing mode theory [T]. 

Previous work describes two different drive mechanisms for the micro-tearing instability. The first is due to the 
thermal force driving parallel currents in the plasma [2]. The second is caused by currents flowing in the collisional 
layers that form around the trapped-passing boundaries, whenever v e /e < lu* [9]. Both of these mechanisms require 
an energy dependent collision operator, and both are dependent on a finite electron temperature gradient. Crucially, 
numerical studies in this paper have shown that neither drive mechanism is particularly important to the MAST 
instability, since the mode is largely unaffected when the energy dependence is removed from the collision operator. 
The trapped particle drive of Catto and Roscnbluth [3] appears to have a small destabilising effect when D e /e < u, but 
there is clearly another drive mechanism for the micro-tearing instability that has not been described in the literature. 

This new drive mechanism requires an electron temperature gradient from which to draw energy, and has some 
dependence on the density gradient. The real frequency of the mode has a roughly linear dependence on both of 
these gradients, which suggests it has drift wave like properties. More detailed studies show the micro-tearing mode is 
strongly destabilised by magnetic drifts and the electrostatic potential. Indeed, if neither of these effects are included, 
then the MAST instability is completely stabilised. 

The MAST micro-tearing mode is unstable in the collisional range 0.05 < 9 e /u> < 2. The lower end of this 
range pertains to the calculations of Catto and Roscnbluth 0] (which have already been mentioned) , while the higher 
end of this range lies somewhere between the collisionless and semi-collisional regimes defined by Drake and Lee 
[7]. Importantly the collisional range of the MAST instability includes the collisional region v e ~ u, which Gladd 
et al and D'Ippolito et al [201 GUI have treated numerically. Nonetheless, despite treating the correct collisional 
range, these papers do not predict the MAST instability since they require an energy dependent collision operator for 
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instability. However, if we remove magnetic drifts from our calculations then energy dependent collisions do become 
more important, thus making some contact with these earlier slab calculations [2T)I |2"T] . 

Importantly, it has been found that the exact form of the ion response has little effect on the MAST instability, 
despite the fact the current layer width d ~ pi. Importantly this does not conflict with previous calculations that 
include the full ion response [8], since these were performed in a different limit (d < pi). 

Finally, the prevalence of micro-tearing instabilities in simulations of the spherical tokamak is not due to its 
characteristic shaping, but rather it seems to be a result of the larger /3, and the smaller radius of curvature in the ST 
(which provides stronger magnetic drifts). 
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Appendix: GS2 Coordinates 

GS2 uses a Clebsch representation for the magnetic field, 

B = VaxVtp (11) 

where ip is the poloidal magnetic flux and a = 4> — q9 + v(9, tp) labels different field lines on a flux surface. Thus, ip 
and a are two coordinates perpendicular to the field in the flux tube, while the distance along the field line is labelled 
using the poloidal angle 9, which extends from 9 = (— tt,tt). In GS2 a and ip are transformed to x and y using, 



aB a p 

dip NORM , 10 n, 
V = ^^p— (13) 

where a is the plasma minor radius, B a is the vacuum magnetic field at the centre of the last closed flux surface, 
ipNORM = ip I '{a? 'B a ), and p n is a normalised flux surface label increasing from zero at magnetic axis to one at the last 
closed flux surface. This results in a coordinate system (x,y,9). 
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